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Abstract. Engineering nano- mechanical quantum systems possessing 
ultra-long motional coherence times allow for applications in ultra-sensitive 
quantum sensing, motional quantum memories and motional interfaces 
between other carriers of quantum information such as photons, quantum 
dots and superconducting systems. To achieve ultra-high motional Q one 
must work hard to remove all forms of motional noise and heating. We 
examine a magneto-nanomechanical quantum system that consists of a 3D 
arrangement of miniature superconducting loops which is stably levitated in 
a static inhomogenous magnetic field. The resulting motional Q is limited 
by the tiny decay of the supercurrent in the loops and may reach up to 
Q ~ 10 10 . We examine the classical and quantum dynamics of the levitating 
superconducting system and prove that it is stably trapped and can achieve 
motional oscillation frequencies of several tens of MHz. By inductively 
coupling this levitating object to a nearby flux qubit we further show that 
by driving the qubit one can cool the motion of the levitated object and in 
the case of resonance, this can cool the vertical motion of the object close 
to its ground state. 
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1. Introduction 

Recently there has been considerable effort towards mapping the boundary between the 
classical and the quantum world by exploring the physics of mesoscopic and macroscopic 
mechanical systems. From an applications point of view, as precision measurement 
of position and acceleration generally involve some kind of motion, the necessity of 
building smaller and more sensitive devices has required a more careful exploration of 
the classical-quantum limit. 

The possibility to couple, control and measure micro-mechanical motion in a wide 
range of different physical systems leads to new experimental applications in different 
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fields such as measuring forces between individual biomolecules [JJ El E], magnetic 
forces from single spins [I] , perturbations due to the mass fluctuations involving single 
atoms and molecules [5], pressure j6] and acceleration |6j, fundamental constants 
[7], small changes in electrical charge [8], gravitational wave detection [9J, as well as 
applications in quantum computation [JU] , quantum optics [TT] and condensed matter 
physics d21[l3]. 

Observing any quantum properties of a mechanical system is a challenge. Under 
typical conditions, energy losses, thermal noise and decoherence processes make it 
impossible to observe any motional quantum effects. To observe quantum mechanical 
motional effects the system has to be close enough to its ground state and it has to 
preserve this quantum coherence for a reasonable amount of time. This leads to the 
necessity of engineering ultra-low dissipative systems (which, in oscillating systems, is 
measured by the quality factor Q representing the energy lost per cycle). To achieve 
this one must engineer a system which is mechanically isolated from it's surroundings 
to an extreme level. On the other hand one must also find a way to cool down the 
motion close to its motional ground state which necessities coupling that system to 
another in order to dump entropy. Numerous nano-mechanical oscillating systems have 
been recently studied such as cavity optomechanical experiments employing cantilevers 
[H], micro- mirrors [To! [To], micro-cavities [TTl ITS], nano-membranes [15] . macroscopic 
mirror modes [20] and optically levitated nanospheres [2Tj (see [22]). As shown in 
Ref. [23], it has been possible to create and control quantum states but, except in a 
few cases, reaching large Q for nano to microscopic sized motional devices is still an 
open problem. In fact, a mechanical oscillator usually involves many degrees of freedom 
coupled together while we are interested in the quantum behaviour of one of them: the 
centre of mass motion. 

In our work, we present a theoretical model for a mesoscopic mechanical oscillator. 
It is inherently non-dissipative (with a estimated Q ~ 10 10 ). It consists of a cluster 
of superconducting loops levitating in vacuum and its motion can be completely 
characterized by the six degrees of freedom of a rigid body. Controlling the system 
is possible through an interaction with a nearby flux qubit which is tuned and driven 
to allow the levitated structure to be cooled close to it's motional ground state. 
We will show how to cool the center of mass motion along one degree of freedom to the 
ground state in an efficient way by coupling the loop inductively to a flux qubit. 

2. Toy Model 

In order to achieve a qualitative understanding of the dynamics of the levitated system 
presented here we will start by studying a simple 2-dimensional toy model [M] which 
describes how it is possible to take advantage of the Meissner effect to obtain a stable 
magnetically levitated mechanical system. For this toy model, we consider a wire loop 
of superconducting material oriented in a vertical plane (Fig. [T| and threaded by a 
discontinuous horizontal magnetic field which is uniform and non-zero above a certain 
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Figure 1. A metallic wire loop in a discontinuous magnetic field. When released to 
fall under gravity the change in flux caused by the loops' motion will induce current 
flow in the wires which will oppose the downwards motion and the loop's height will 
undergo harmonic oscillation. For wires made from normal metals the oscillations 
are damped due to resistance within the wires and the loop will eventually fall but 
for superconducting wires the oscillations will continue for many periods yielding a 
motional oscillator with extremely high motional Q factor. 



height and which vanishes below. Initially the loop is at rest at a position where it 
crosses the magnetic field discontinuity line. 

When the loop falls under the force of gravity, the changing magnetic flux within 
the loop will induce circulating currents which generate Lorentz forces that counteract 
the downwards motion of the loop. If the loop is made up of normal conducting material 
the loop will undergo harmonic motion in the vertical plane which will be damped (due 
to resistances within the wire). The currents generating the Lorentz forces will dissipate 
and the loop will eventually fall down. However, let's consider the wires to be fashioned 
from superconducting material. In this case things change due to the Meissner effect 
which induces currents such that the flux <p through the loop remains constant at all 
times. With reference to the notation set up in Fig. [TJ the flux at loop position x 
is given by 4>(x) = wBx. The flux due to the induced current I(x) is given by the 
usual formula LI, where L is the loop inductance. Thus, the Meissner effect tells us 
that: LI(x) = wBx. Using this result and noting that the loop's potential energy is 
V(x) = \LI 2 (x) — mgx, one can see that this system corresponds to a gravitationally 
driven harmonic oscillator. 



3. Model 



The above 2D toy model is obviously not completely trapped if naively generalised to 
3D. To achieve stable three dimensional trapping we must generalise the above toy 
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model in a more sophisticated manner. Let us now consider the system sketched in 
Fig. [2] We consider a spatially strongly inhomogeneous magnetic field generated by a 
fixed magnetized sphere placed above a cluster (the motional resonator) of three non- 
intersecting loops oriented along mutually orthogonal planes. This particular shape will 
allow the potential energy of the system to have a stable minimum and also leads to 
a diagonal inductance matrix, so that the current flowing in one loop does not induce 
any currents in the other loops. A flux qubit is placed at a certain distance from the 
center of the cluster of loops. The inductive coupling between the flux qubit and the 
cluster will provide us with the control required to establish motional cooling of the 
levitated cluster. In the following we first study the system classically to establish full 
3D trapping, translation/rotation oscillation frequencies at the trap minimum, before 
studying the quantum mechanics and motional cooling protocols to cool the vertical 
translational degree of freedom towards the ground state. 

The total Hamiltonian H is the sum of terms involving the motional resonator (H r ), 
the qubit (H q ) and the interaction (Hi): 

H = H r + H q + Hi , (1) 

and in the following sections we develop each of these individual Hamiltonians. 

3.1. Resonator 

The Hamiltonian of the cluster of loops (which we will now denote as the motional 
resonator or simply resonator), can be written in the following form: 

H r = K CM + K rot + V , (2) 

where Kqm is the kinetic energy due to the translational motion of the center of 
mass, K rot is the rotational kinetic energy and V is the effective potential energy as a 
function of the translational and rotational degrees of freedom. We will now define some 
notation used in the following. 

• m: mass of the resonator. 

• X: vector position of the various part of the resonator in a co-rotating reference 
frame with origin at the center of mass. 

• R: coordinates of the center of mass in the laboratory reference frame. 

• r: vector position of the various part of the resonator in the laboratory frame. 

• r = f — R 

Note that the origin of the laboratory frame is taken to be at the center of the magnetic 
sphere (see Fig. [2]) . 

Because of the rigid body properties, the only allowed relation between r(X,t) and 
X is given by r(X,t) = 0(t)X, where O = e^* ai ® Ti (with (Ti) jk = e^) is a ro- 
tation matrix and oli G [— 7r,7r]. We can then define the angular velocity vectors as: 
Qi = 0(t) _1 0(/f:) = d/Tj. The motion of the rigid body is thus completely determined 
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Figure 2. Cluster of three insulated superconducting loops levitating in a magnetic 
field generated by a magnetized sphere with magnetization vector directed towards the 
positive i direction. (A) In the laboratory frame the axes are labeled {xj}j =1 . For each 
axes we have angular coordinates a L (labeled accordingly) corresponding to rotation 
around that axes in a counterclockwise sense. The reference system has the center 
of the sphere as origin. The loops are considered non-intersecting (i.e. no electrical 
conduction between the loops). The dimensions of the sides of the cluster of loops 
parallel to axis Xi we take to be hi. We depict the magnetic vector field generated 
by the spherical magnet and the nearby flux qubit (which we take to be a 3-junction 
phase qubit), sitting under the cluster on a yellow substrate. (B) More detailed view 
of the cluster of loops. 



by the set of variables (R,a). 

Given these definitions, the inertia tensor of the resonator is ly = f dV r p(X) (^X 2 5ij — Xi ■ Xj 
The kinetic energies are defined with respect to the reference system in Fig. [2] as: 
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The potential V is just the sum of the flux energy due to the current flowing in the 
loops and the gravitational potential energy: 

1 3 

V=-Y J L a I 2 a -mgR l , (5) 

0=1 

where the index a = 1,2,3 labels the loops accordingly to the direction 
perpendicular to the plane of the loop (where 1 — > i, 2 — >• j : 3 — > k), L a is the inductance 
of the loop a and I a is the current flowing in the loop a. By the symmetry of the resonator 
the mutual inductances between the loops is zero. The currents are just functions of the 
six degrees of freedom (translations, rotations), describing the motion of the resonator 
via equations originating from the Meissner effect: 

A(f> a (R,a) + LJ a (R,a) = , (6) 

where A(f> a (R, a) = 4> a (R, a) — 0a(-R(O), a(0)) is the difference in magnetic flux threading 
loop a when the system is in the configuration labeled by (R, a) and when the system 
is in its initial configuration (-R(O), c?(0)). Any infinitesimal change in flux due to an 
infinitesimal displacement/rotation induces a supercurrent whose action is to restore the 
loop's position/orientation. The stronger this restoring force is the higher the oscillation 
frequency will be. We now compute the dependence of the magnetic flux on the system's 
configuration. 

We start by calculating the magnetic flux threading the loops. The vector potential 
A generated by a sphere with homogeneous magnetization vector M calculated at the 
point r + R (where R is the center of mass position vector referenced from the centre 
of our coordinate system - the centre of the sphere), is just: 

A(f ] M,R) = ^-±^-MA(r + R) . (7) 

If we denote with E a the area vector of loop a, then the flux through this loop can 
be expressed as: 

4>a(R,oi)= [ B-dt a = [ A-dr . (8) 

In a first order approximation and by supposing that the initial position of the 
resonator is such that a, = 0, we have: 

A0 a = ^{K a ■ [f- R(0) A a] + Q a ■ [M A a}) , (9) 

where the vectors K a and Q a are easily calculated from the magnetic field and the sphere 
magnetization. By using Eq. Q and Eq. @ we can expand the potential energy to 
second order in spatial and angular coordinates as: 

E = l -E i3 CC 3 , (10) 



8 



where ( = {x,y, z, a x , a y , a z }, is a vector whose components are all the parameters 
describing the motion and where the matrix elements depend on the vectors defined 
just above. 



In Section [5] we will show how certain choices of the parameters of our system can 
lead to a decoupled degree of freedom (the vertical direction with associated variable 
Ri), with a quite high oscillating frequency u r . This frequency is also well resolved 
with respect to the frequencies of the other modes (the difference between frequencies 
is bigger than the linewidth of the resonator) , hence it can be well resolved by the qubit 
coupling. Note however, that this mode is strictly decoupled from other rotational and 
motional modes only to second order in perturbation theory hence it is important that 
the system be initialized not too far from equilibrium. We will consider the quantized 
variable corresponding to small deviations from the initial position along the i axis: 
x = Ri — x (x being the x component of the vector R ), whose quantum operator we 
define as: 



r l m — — )(a + a t ) . (11) 
2muj r 

The dynamics of this variable is governed by an harmonic Hamiltonian whose frequency 



co r can be obtained through the energy given above in Eq. (10). 



3.2. Qubit 

The flux qubit is considered to have characteristic frequency u q and to be driven by 
a classical field of frequency u d slightly detuned by u q . The detuning parameter is 
5 = Ud — u q . The field drives the qubit with Rabi frequency Q. The Hamiltonian in the 
rotating frame with frequency Ud is: 

H q = --a z + -a x . (12) 

3.3. Interaction 

The classical interaction Hamiltonian for the inductive coupling between the loops of 
the cluster and the loop of the flux qubit is given by: 

Hi = M rq I r I q , (13) 

where M rq is the mutual inductance between the qubit and the loop of the resonator 
which is horizontal in the initial position and I r = I\ and I q are the currents flowing 
in that loop and in the qubit. Because we are considering small deviations from the 
equilibrium initial position, we neglect the mutual interactions between the other loops 
or due to the angular motion of the object. We will consider the dependence of the 
coupling in the vertical direction by performing a first order expansion. We expand I r 
to first order in small deviations from the initial position using d9l). In the quantized 
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version we also replace /„ with I„a z (see for example [213 )• Denoting Dj(0) = ^ as 

ax R{0) 

the derivative of the current evaluated at the initial cluster position and remembering 



(11) we write the quantized interaction Hamiltonian as: 



Hi = h ^{a + tf)& z , (14) 

with: 

The mutual inductance is simply calculated as: 

M -^J dJ f' < 16 > 

where ds and ds 1 are vectors tangent to the two loops and R is the distance between 
the infinitesimal part of each loop and the integration is along each of the two loops. 

The final quantized expression for the Hamiltonian describing the motion of the 
cluster in the vertical direction and its interaction with the qubit is (setting h = 1): 

S Q A 
H = --a z + —a x + u r d ] d+ -(d + d ] )a z . (17) 

In order to be in the motional quantum regime we must devise cooling schemes to bring 
the resonator close to its motional ground state as described in the next section. 

4. Cooling 

4-1. Cooling scheme A: Introduction 

Following [26] we will consider the resonator and the qubit coupled to separate thermal 
baths and interacting with each through the coupling Hj (see Fig. [3]). 

The system's quantum motional state, p evolves under the following master 
equation: 

j5 = -i[HJ)+C r (p) + C 7 (p) . (18) 
The qubit Liouville operator contains decay and dephasing with rates T± and IV 

£r(p) = ~Y^ N q + l)(2£-po+ - or + a_p - po + d_) + 
+ -^( Ar <?)( 2 ^+P^- - v-v+P ~ P<5--<7+) + 

+ ^(^pa z -p) , (19) 
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where N q = (e k B T — 1) 1 ! and where T is the temperature of the qubit bath. The 
Lioville operator for the motional resonator is: 

£ 7 (p) = —(N t h + l)(2apa* — a 1 dp — pa 1 a) + 

+ ^{N th )(2d ] pa - dd ] p - pdd ] ) , (20) 
where the mechanical dissipation factor 7 has been introduced, and the equilibrium 

hu> r 

phonon occupation number is Nth — (e fc s T — l) -1 . 

We now trace out the qubit to find an effective equation of motion for the resonator. 
From this we will see that under certain cases we can obtain a motional cooling process 
bringing the cluster towards its motional ground state. We will consider the two different 
temperature regimes studied in [26] . 




Figure 3. Motional Cooling Scheme: we inductively couple (with strength A), the 
motion of the levitated resonator with a nearby flux qubit. The resonator looses 
energy at rate 7 due to coupling to a local phonon bath while the qubit experiences 
both decay and dephasing at rates T± and Tu. In addition we drive transitions in the 
qubit at rate O. 

4-1.1. Low Temperature Regime Here we consider a regime where A a/ N t h + 1/2 <^ 
r±,u r , which is equivalent to the Lamb-Dicke regime in laser cooling (see [26] and ref- 
erences within). The first condition tells us that the qubit damps much faster than 
energy can be transferred to it from the resonator. By observing that the parameter A 
describes the position displacement kick of the resonator when the qubit undergoes a 
spin flip (see [22]), the second approximation essentially says that the back-reaction on 
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the resonator following a qubit flip is small relative to the oscillator energy. 

Thanks to these assumptions one can write the full density matrix as p(t) ~ 
Pg <S> p r (t) where p° q is the qubit steady state density matrix. In this way one can obtain 
an effective master equation for the resonator after having traced out the qubit. The 
result can be encoded in a new effective mechanical damping factor given by V = Tc + 7 
with r<7 = S(uj r ) — S(—uj r ) where S(u) denotes the qubit fluctuation spectrum: 

\ 2 roo 

S(u) = -Re J e^dr {(a z (r)a z (0)) ~ , (21) 

and where (-) denotes the steady state expectation [26J. What interests us is the 
final motional occupation number which is given by: 

n LD = ^ + N , (22) 
t c 

where iV = S(— u r )/Yc, and the subscript stresses the Lamb-Dicke regime we are 
supposing. 



4-1.2. High Temperature Regime As one can see from (22), for low enough N t h the final 
occupation number is given by N . Increasing the temperature of the motional bath we 
see that the final occupation number grows linearly with Nth- However, when Nth is 
so large the previous assumptions A a/ N t h + 1/2 <C T±,u r do not hold true anymore it 
turns out [26] that, under the temperature independent assumptions (A T±,co r ), one 
can study situations in which the final occupation number is given by the more general 
formula: 

n f = N ^ (^+ i + C eV(W) ) > ( 23 ) 



with parameters defined in Appendix A As described in that section, in order to obtain 
this formula one must consider the resonator's initial quantum state is given by the 
coherent state with \a\ 2 ~ Nth a coherent state labeled by the parameter a as 

initial condition. Given this assumption the "renormalized" cooling rate Tc is an a 
dependent quantity. 



4-2. Cooling scheme B 

Rather than use the scheme depicted in Fig [3j one can also cool the motion of the 
resonator (at least in the case of near resonance with the qubit) by performing a sequence 
of repeated measurements on the flux qubit inductively coupled to the resonator's 
position [27] . 

Using a flux qubit made out of three superconducting loops with four Josephson 
junctions one [2HI EH] can tune the magnetic energy bias HS to zero. We also assume 
that the inductive coupling strength is much smaller than the qubit frequency (A u r ). 
This will allow us to perform a rotating wave approximation on the Hamiltonian of the 
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full system. We consider the system to be initially prepared in the separable state 
Po — \g)(g\ ® p T with the resonator in the thermal state = (n\p r \n) = ■ Then, 

one performs N measurements on the qubit (with associated operator given by \g)(g\ 
and \e)(e\ where g (e) is the ground (excited) state) at times tj = jr where r is short 
time interval. One finds that, provided the measurement outcome is always |<7)(<7|, at 
the conclusion of these N measurements the state of the resonator is [27] : 

MAO-E ^™'" 1 . (24) 

n>0 ^ ' 

where /i = e iAr/2 , /i n = e- i(r ^ 1/2) ^ T (cos Vt n r + 2sinfi n rcos2# n ), tt n = 
^ (A — cov) 2 /4 + g 2 n (n > 1), tan2# n = 2gy / n/(A — uj n ) (n > 1) and where N indicates 
the number of measurements. P{N) is the probability this sequence of measurement 
outcomes occurs. 

5. Results 

Here we consider our 3D model in a more realistic setting and adopt physical parameters 
in the table shown in section [8} Using this physically realistic set of parameters we now 
examine in more concrete detail the classical and quantum mechanics of the resonator 
and the levels of cooling achieved by the above outlined methods. 



5.1. Classical Motion 

We now examine the classical motion of the resonator by focusing in particular on its 
stability and its frequency of oscillation. 

The effect of the gravity in the potential energy V(x) causes a small and negligible shift 
of the potential minimum. 

By considering the magnetization of the sphere to be aligned along the vertical x direc- 
tion and expanding to second order the potential energy around its minimum, one finds 
that the x-motion decouples from the other directions. One finds further the existence 
of three degrees of freedom which are "zero modes" , i.e. do not contribute to the energy 
up to third order. One zero-mode is related to the rotational symmetry around the 
vertical axes. The other two are due to the coupling between the variables y(z) and 
a z (a y ). however, if one considers these contributions to higher orders these degrees of 
freedom contribute to the energy value to give overall stability as shown in Fig. |5j For a 
cluster of loops of "typical "dimensions (as defined below) of (1, 10, 10) \xm and with a 
wire thickness of 0.1 |^m the root mean square thermal motion of the cluster at 15 mK 
is ~ 10 -12 m along the spatial directions and 10~ 6 radians along the a y and a z angular 
directions (while there is symmetry along the other angular direction). 

The frequency of the translational motions depends on the physical dimensions of 
the cluster of loops. In Fig. [6] the frequency in the x and ^-direction is plotted as a 
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Figure 4. Exact numerical evaluation of the potential energy (in units of Eq = hui r ) 
as a function of the vertical translational degree of freedom x (with origin the stable 



point and in unit of xq — 
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Figure 5. Exact numerical evaluation of the potential energy as a function of the 
horizontal translational degree of freedom y (with origin the stable point) and the 
angular degree of freedom a z . Because the potential is confining to third order in 
these degrees of freedom the overall system is stable. 



function of a scaling factor k for the system. The values for all the parameters are 
explicitly shown in the table in section |8j 

The frequencies decrease as the scaling grows, showing that the trapping in the transla- 
tional degrees of freedom is becoming less tight, but still present for a large range of size 
scales. In regards of the rotational degrees of freedom, even when the system is scaled 
up with a factor 10 with respect to the typical parameters, numerically we find that the 



trapping is still achieved in the small oscillation regime (see section 3.1). 



In summary, using the parameters shown the table in section [8j for small loops 
and very high magnetic field inhomogeneities the cluster of loops experiences stable, 
levitated magnetic trapping. The frequency of oscillation along the x-direction is larger 
than in the other directions making it easier to cool. For these reasons, in the following 
we quantize the x-translational motion. 



14 




20 40 60 80 100 



Scaling Factor (K) 

Figure 6. Translational mode frequency as a function of system size. The frequency 
along the x direction (blue line) and y (and for symmetry reasons z) direction (red 
line) are plotted as a function of the scaling factor k which defines the dimension of 
the system. The radius of the magnetized sphere is defined as R sp — k \xm. The 
dimension of the parallelepiped enveloping the cluster of loops is k x (0.1, 1, l)(j.m and 
the thickness of the wires is k x 0.01 urn. The cluster mass scales as k 3 . When 
considering different values for the scale factor, the initial position of the center of 
the cluster of loops is chosen to keep the distance between the center of the sphere to 
the top of the cluster fixed at ~ 1.1 R sp where R sp is the radius of the sphere. The 
frequencies are higher for smaller sizes of the system. As "typical" parameters for the 
system we will consider the ones given by a scaling factor k = 10. 

5.2. Superconductivity and Dissipation 

The superconducting loops are considered as made out of Niobium-Tin (NbSn) which, in 
normal conditions, has very high critical magnetic field strength (~ 30 T ), although, in 
our case, the situation is more critical due to the small thickness of the superconducting 
wires. 

We now examine two types of motional decay mechanisms and estimate the resulting 
motional Q for the resonator's motion. We first consider the loss in energy from the 
moving resonator due to its action as a dipole emitter of radiation. 
In Fig. [7] the intensity of the magnetic field (in Tesla) is plotted as a function of distance 
from the center of the magnetized sphere (radius 10 p. m). 

As we saw above, the motion of the cluster of loops is trapped in all three directions. 
Due to the Messnier effect, the greater the spatial excursion from its equilibrium position 
the greater the current density inside the wires. Considering the physical parameters 
given in the table in section|8]we find typical current densities of ~ (10 8 , 10 8 , 10 7 )A/m 2 
(respectively in the loops perpendicular to the x, y, z directions). Considering each 
current to be described by I = Ioe~ tult , where u is the top frequency in each direction, 
by approximating the loops to be circular we can calculate the power dissipated by the 
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Figure 7. Magnetic field intensity along the vertical x. 



oscillating currents as electromagnetic radiation as (see for example [30J ) : 

1 n An lixr Iq 




(25) 



where r is the radius of the loop and A = cjv (where c is the speed of light). From 
this we observe that the power radiated is absolutely negligible. 

However, the main source of dissipation is the viscous drag of flux lines oscillating inside 
the pinning wells inside the superconducting wires [7]. One can estimate the motional 
quality factor Q = — for a system similar to (the two dimensional version of) the one 
presented here to be Q ~ 10 11 . In our calculations we will consider Q 10 10 . Such 
an enormously high motional quality factor is already known in levitating systems (31] 
(Q ~ 10 12 ), although typically one finds motional Q values in the range ~ 10 3 /10 6 for 
cavity optomechanical experiments [22]. The prospect for ultra-large motional Q is one 
of the primary benefits of our scheme. 

5.3. Coupling between the resonator and the flux qubit. 

In this section we will study the coupling between the resonator and the flux qubit. The 
possibility to tune this coupling arises from the (distance dependent) inductive nature 
of the coupling. 

Fig. [8] shows the value of the coupling strength between the qubit and the resonator as 
a function of the distance d between the center of mass of the resonator and the center 
of the loop of the flux qubit the latter taken to be a circle of radius 5|xm. 

We observe that the coupling strength can be adjusted over quite a large range by 
fixing the distance between the qubit and the resonator. 

5-4- Cooling of the motional resonator. 

5.4-1. Cooling scheme A. In this cooling scheme the motional energy of the resonator 
is damped away thanks to the faster decay of the flux qubit. 

We will start by studying the low temperature limit. This means that the N t h of the 
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Figure 8. Coupling constant between the qubit and the resonator. The coupling arise 
from the mutual inductance between the loops of the resonator and the qubit loop. 
By supposing the loop of the qubit to be placed in a plance orthogonal to the vertical 
direction, the contribution to the coupling comes, in a good approximation, from the 
loop of the resonator perpendicular to the x direction. The distance d is defined as the 
distance between the center of mass of the resonator and the center of the qubit loop. 
The qubit loop is supposed to have a 5p.m long radius. 



bath connected to the resonator must be such that X^J N th + | <C T±, oj r is fulfilled. 
In Fig. [9] we plot the final motional average Fock number as a function of the qubit 
resonance frequency u q and Rabi frequency Q, while in Fig. 10 we plot the same final 
occupation number as a function of u q after fixing Q to an optimal value. In particular, 
we chose the Rabi frequency f2 and the detuning 5 to fulfill the relation uj q = \/Vt 2 + 8 2 . 
For comparison, the dashed line in Fig. [9] shows the equilibrium phonon number if the 
resonator is in equilibrium with a 15 mK thermal bath. The cooling happens in the 
region defined by the dashed line. Fig. 10 , which is a section of Fig. [9] obtained by fixing 

n. 

Now we move to a regime where no constraints on the temperature of the bath 
connected to the resonator are imposed so that the only relations to be fulfilled are 
A <C r_|_,a; r . This will allow us to use an higher coupling constant with respect to the 
low (but still non-zero) temperature limit, previously considered. 

Fig. \TT\ shows the final expected cooled final motional Fock number n r in the case where 
there exist a large frequency mismatch between the qubit and the resonator. As expected 
the cooling is quite poor. For comparison, we also plot the dashed line labeled as tild 
obtained from the low temperature theory extrapolated to high temperature (as given 
by Eq. (22)). One can see that in the low temperature regime the two theories agree 
while at high temperatures they disagree. Above a certain bath temperature cooling is 
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Figure 9. Contour plot representing the final occupation number (blue corrsponding 
to zero phonon and red corresponding to 200 phonons) versus the qubit parameters fl 
and Lj q with 5 fixed by the constraint S = ^/cj^ — Q 2 . The dashed line represents the 
phonon equilibrium number if the system is in equilibrium with a thermal bath with 
temperature 15 mK. Cooling takes place for parameters inside the dashed line. The 
frequency of the resonator has been taken to be such that oj r = 27rl0 7 . The coupling 
between resonator and qubit has been set to 10 4 . With this value, the low temperature 
approximation is a good one for qubit frequencies such that uj q > 10 9 . 



no longer possible. Below this bath temperature cooling is possible though in the far 
off-resonant case the degree of cooling is poor. 



In Fig 11 we consider the case where the qubit and the resonator are vastly off resonance 



and this gives very poor cooling performance. One can expect the cooling to be far better 
when they are both on resonance. To achieve this, one must either reduce the qubit 
frequency or increase the resonator frequency, possibly via new flux qubit designs [32] or 
by replacing the magnetized sphere with a magnetic tip [33J with much larger magnetic 
field gradients. 



In Fig. 12 we consider the same cooling process as in Fig. 12 but with the resonator 
now resonant with the qubit. In this case the system shows a significant enhancement 
in the cooling rate which scales from 7 = 0.1 (off-resonance) to 7 + Tc G (5 x 10 2 , 10 2 ) 
(resonant) for the regime where the cooling is still possible. The net result is that, even 
if the resonantor is connected to a very hot thermal bath (Nth ~ 10 4 ), one can cool the 
system to a final average phonon number around rif ~ 1. 

We showed that the interaction between the qubit and the resonator leads to 
an increase of the cooling rate. Cooling is possible even when the initial effective 
temperature is quite high. Cooling efficiency increases significantly when the resonator 
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Figure 10. Plot of the final occupation number as a function of ui q with qubit 
parameters fixed as follows: il — 0.950;^ and 5 = ^Juj^-Q? . This graph is just 
the cross section of fig. [9] obtained by fixing J7. The dashed line is the same as in fig. 
[9] and represents the phonon equilibrium number if the system is in equilibrium with a 
thermal bath with temperature 15 mK. The frequency of the resonator has been taken 
to be such that u> r — 27rl0 7 . The coupling between resonator and qubit has been set 
to 10 4 . With this value, the low temperature approximation is a good one for qubit 
frequencies such that u q > 10 9 . 

10 6 
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JVth 

Figure 11. Plot of the final occupation number as a function of the initial occupation 
number. The dashed line labeled as uld represents the low temperature limit. The 
other dashed line represents the identity line. For this plot ui r — 27rl0 7 , u> q = 10 9 , 
A = 10 5 . 

and the qubit are in resonance. 

5.4-2. Cooling scheme B. This scheme cools down the motion of the qubit using 
repetited measurements on the qubit inductively coupled to the resonator. The cooling 
is achieved in the case in which the outcomes of the measurements always tell us that 
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Figure 12. Plot of the final occupation number as a function of the initial occupation 
number. The dashed line labeled as uld represents the low temperature limit. The 



other dashed line represents the identity line. Here lo t = 10 



10 9 , and A = 10 5 



The result is that cooling is possible in this resonance situation even if the bath 
connected to the resonator is very hot (T ~ 100 K). 



the qubit is in the ground state. This "lucky" sequence happens with probability 
that decays exponentially with the number of measurement but which, in the limit 
of infinite measurement does not vanish, but rather converges to the initial ground state 
probability. 

Henceforth we suppose the qubit and the resonator are in resonance. The value of the 
frequency has been fixed so to maximize the results given a maximum value for the 
coupling A = 5 MHz. Fig. 13 shows the final occupation number as a function of the 
number of measurements on the qubit. In this way it is possible to reach the ground 
state after a quite limited number of measurements. As explained in [27], requiring equal 
times between measurements (which would be a technical difficulty) is not important 
and can be relaxed by introducing random time deviations. The situation would even 
improve in this case. 



6. Conclusions 

We presented here a mechanical oscillating system that is inherently non-dissipative 
that utilizes superconducting material levitated in a vacuum via the Meissner effect. 
We showed that the cluster of loops is trapped in a potential well given by a magnetic 
field generated by a nearby magnetized sphere. Among the total of 6 degrees of freedom 
of the resonator we quantized the variable associated with the vertical direction as, 
in a second order approximation, it is decoupled from the other degrees of freedom 
and it possess the largest oscillating frequency. We also showed how to cool down the 
motion of the system in the low temperature regime as well as in the high temperature 
regime through coupling to a qubit. The system allowed us to consider two different 
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Figure 13. Plot about the cooling by mean of measurements on the qubit. The blue 
line represents the plot of the final occupation number fij as a function of the number 
of measurements on the qubit. The red line represents the survival probabilty P(N) 
which, for N — > oo, goes to the initial probability to find the system in the ground 
state. For this plot the resonator is interacting with a T = 15mK thermal bath and 
we consider a resonant situation with u r = u q = 2.5 x 10 8 Hz, A = 5 x 10 6 Hz and 

T = & S. 



cooling schemes which work very well in the case resonator and qubit are resonant. The 
system can be improved by considering stronger or higher gradient magnetic fields, which 
could be achieved by replacing the magnetized sphere with other geometric magnetized 
objects. This is not the only direction for further studies. In fact a complete analysis 
would quantize all the degrees of freedom and would find a way to cool all of them (for 
example by introducing other inductive couplings with other flux qubit). 
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8. Parameter Table 



Variable 


Reference Value 


Meaning 


R 


10 |xm 


Radius of the Sphere 


Ro 


110% R+lhi 


Initial distance from center of sphere 


M 


— m~ 3 

Ho 


Magnetization density 


M 


X 


Magnetization vector 


P 


6.6 x 10 3 kg m~ 3 


Superconductor mass density (close to Nb 3 Sn density) 


r 


50 nm 


Wires radius 


hi 


l|j.m 


Wires length along x 


h 2 


10|J.m 


Wires length along y 


h 3 


10(J.m 


Wires length along z 


d 


> l^m 


Distance between qubit and cluster 


n 


(0.3 - 10) GHz 


Driving Field Frequency 


5 


-(1 - 100)%Q 


Dephasing 


Ti 


0.5m.s 


Qubit Decay Time 


T 2 


l\xs 


Qubit Dephasing Time 


1 _L 


1 1 

Ti 2A r T LS+i 


Qubit Decay Rate 




T 2 l - |Tf 1 


Qubit Dephasing Time 


A 


(10 4 , 10 5 )Hz 


Interaction between qubit and cluster 


T 


15 mK 


Temperature 


Q 


10 10 


Estimated Quality Factor 




2tt10 MHz 


Mechanical Angular Velocity 


7 


n Q 


Mechanical Energy-Damping Rate 






Thermal Equilibrium Occupation Number 


% 


10 nA 


Qubit current 


m 


4 x 10" 15 kg 


Cluster of loop mass 




w 7 x 10 5 A/m 


Current derivative calculated at Origin 



Appendix A. High temperature regime final occupation number 

By supposing the initial state of the resonator to be a coherent state parametrized by 
the complex variable a, and by supposing the time scale of the interaction to be much 
longer than the time scale with which the resonator interacts with its bath (ATi 1) 
one can obtain, from the master equation of the system, the bloch equation given by: 

< S >= A < S > -i\[e~ iwA a + e iu>rt a*)A <S> -T ± V Z , (A.l) 
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— * 

where S = {<r_, <r +i , a z } represents the pauli operators acting on the qubit space and 
where A and V are given by: 



/ i8-± 



\ 





in 



-is - 1 




(A.2) 



and: 



V z = 





V 1 / 



By writing the solution in the following representation: 



+oo 



<S>(t)= Sn(t)e-" 



— inu) r t 



(A.3) 



and by supposing that S n {t) = S n (no time dependence), one can write the steady state 
bloch equation in the following way: 



with: 



and: 



T^n-i + B n S n — T*S n +i — V n , 



T = a\A , 



B n = A + iu r n , 

— * — * 

V n = S n F±V Z . 



(A.4) 

(A.5) 

(A.6) 
(A.7) 



To this equation introduce a parameter m G N and solve iteratively in m supposing, 
at each step, that Si — for % > m. For what follows we are interested in the solution 
for the terms Sf and Si ± . They are given by the following continuous fractions: 



-T*S , 



K i 

S+i = -r r — TS'o , 



K 



+i 



(A.8) 
(A.9) 



with: 



(A.10) 
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K i 

K+i 
R 



■ 1-2 
1 



-T* - T + B +1 

•*l +2 

K_i + £ + K+i , 



(A.ll) 

(A.12) 
(A.13) 



where • and •* are denned recursively as: 



• 1-3 
1 



T* + B. 



In 



-T* - T + B 
■ * i 

1+3 



+2 



(A. 14) 
(A.15) 



One can stop after a certain amount of iteration by simply substituting, instead of 
the previous recursive relation, the ending one: 



•|-n = B -n , 



In 



B 4 



By following [26] the final occupation number is given by: 



/ 



rif = N th 



e N thV 2 



J 



V 1+ c 

where the cooling rate T{a) is given by: 



(A.16) 
(A.17) 



(A.18) 



T(a) = iX 



a or 



with Sf and are given by the solutions A. 8 and the other terms are given by: 



(A.19) 



C 
h 



r c (Q) 

7 

da aY r ( r 

r c (a) 

r c (o) ' 

A 



(A.20) 
(A.21) 
(A.22) 
(A.23) 
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